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We consider general-relativistic rotational perturbations in homogeneous and isotropic Friedmann 
- Robertson - Walker (FRW) cosmologies. Taking linear perturbations of FRW models, the general 
solution of the field equations contains tensorial, vectorial and scalar functions. The vectorial 
terms are in connection with rotations in the given model and due to the Sachs - Wolfe effect 
they produce contributions to the temperature fluctuations of the cosmic microwave background 
radiation (CMBR). In present paper we obtain the analytic time dependence of these contributions 
in a spatially flat, FRW model with pressureless ideal fluid, in the presence and the absence of a 
cosmological constant A. We find that the solution can be separated into an integrable and a non- 
integrable part as is the situation in the case of scalar perturbations. Analyzing the solutions and 
using the results of present observations we estimate the order of magnitude of the angular velocity 
corresponding to the rotation tensor at the time of decoupling and today. 

PACS numbers: 98.80.Jk, 98.70.Vc, 98.65.-r 



I. INTRODUCTION 

The importance of the study of possible contributors to the temperature fluctuations of the cosmic microwave 
background radiation is indisputable. The large angular-scale anisotropy is dominated by the Sachs - Wolfe effect 
where the first order perturbations in the metric tensor give rise to the appearance of temperature fluctuations in the 
CMBR. 

In present paper we focus on the effects of rotational perturbations in a flat FRW model filled with a single, 
incoherent, perfect fluid matter source, and we analyze its contributions to the fluctuations both in the absence and 
the presence of a cosmological constant. This problem certainly has been studied by several authors before, eg. 
0; S 0| ) however none of these works payed too much attention on the possible analytic solutions in the presence of 
A. Throughout this paper we follow a systematic treatment in calculating parallely with vanishing and nonvanishing 
cosmological constant and comparing the obtained results of the calculations. In Section II. we take a general overlook 
on the C°° solutions of the perturbed field equations given by Sachs & Wolfe flj. White [5] and Perj es et al. m. and the 
rotation tensor. In the presence of the cosmological constant we use the physical time coordinate t in the calculations 
- as the background solution has a closed analytic form only with this time variable - while in the absence of A we 
use the conformal time coordinate 77. In Section III. we calculate the contributions of rotational perturbations to the 
temperature fluctuations due to the Sachs - Wolfe effect and analyze the obtained expression. We also give a recipe 
in this section, how to calculate the anisotropy power spectra of the temperature fluctuations caused by rotational 
perturbations. In Section IV., using the result of present observations, we give an upper estimation for the order of 
magnitude of the rotations and the corresponding angular velocity at the time of decoupling and today. 



'Electronic address: 'czinner@rmki.kfki.huj 
tElectronic address: vasuth@rmki.kfki.hu 



2 



II. LINEAR PERTURBATIONS 



In the case of a perfect fluid matter source the Einstein equations are 

Gl = {p + p)u'^u,-p5l , (1) 

with the four velocity u° = dx'^/dt of the fluid normalized by u'^Ua — 1, p is the pressure, p is the energy density and 
is the Einstein tensor. The line element in a spatially flat FRW spacetime with Cartesian coordinate system is 

ds^ gabdx"'dx'' = a'^ {r])r]ahdx'' dx'' , (2) 

where rj is the conformal time coordinate with the definition dt = adrj, a{ri) is the scale factor and rjat is the Minkowski 
metric with the following signature, 

Vab = rf' = dtag{+l, -1,-1,-1). (3) 

Here and throughout this paper Roman indices run from to 3, Greek indices from I to 3 and we employ units such 
that c = SttG = 1. 

The well-known unperturbed solution of the Einstein equations in the case of incoherent matter without cosmological 
constant is 

Hr 7]^ 3Hr 

where Hr is the present value of the Hubble parameter H = a! ja^ , the prime denotes the derivative with respect to 
the conformal time and we chose the r\R — 1 normalization condition. The also well-known background solution of 
the field equations for incoherent matter in the case of A ^ is 

p = , a = aosinh2/3(Ci + C^) , pa"" ^Cm . C , v= f- , ^fl = 1 , Co = , (5) 

2 J a 

where Qq is restricted by the normalization condition of rjR. 

Hereafter in this section we follow the perturbation method of Sachs & Wolfe [l| and White Q and we write the 
perturbed metric gab as 

ds^ = gabdx^-dx^ = a^(r]){f]ab + hab)dx''dx^ , (6) 

where hab is the metric perturbation. We use the Minkowski metric i^ab and its inverse to raise and lower indices 
of hab and of other small quantities. In the perturbed spacetime we choose comoving coordinates, where the four 
velocity of the fluid and the metric perturbation /iqo are 

= ^ , /loo = . (7) 

a 

The general C°° solution of the perturbed field equations for incoherent matter without cosmological constant is 
given by Sachs & Wolfe Q and White [|, 

ha/S = -TT + ^ ? l^Q,/3 + ^/3,a) ^ ^ + Valii^ - —^ti ap , (8) 

rj orj \ rj J \ V V J V 10 

hoc. = -^V^Ca , (9) 

where and throughout this paper comma denotes partial derivative and A = A(.t^), B = B{x^) are two scalar 
functions corresponding to the potentials for density perturbations. The Ca = Ca{x^) vector function - satisfying 
the C"q, = transversality condition - is related to the perturbed rotation (vorticity) tensor LOab, see (HZ]). The 
Dap{x^ ,7f) transverse, trace- free tensor function represents gravitational waves as the solution of the following wave 
equation: 

^ - V^^ D^p^Q, D^p = Dp^ , D^p^ = , Z?'; - , (11) 
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where = A is the standard Laplacian defined by V^/ = —r^^^^f^^^, for an arbitrary smooth function /. 

The general C°° solution of the perturbed field equations for incoherent matter in the case of A 7^ is given by 
Perjes et al. 0, 

Kp - S^f, + {^^^^ - coth(Ct)^ {Cc.,p + C/3,a) + A,„/3 coth(Ci) + T^o^pB - i^^i^^^S,^^ , (12) 

= ~sinh^/\ci)^'^" ' ^''^ 
Cm cosh(Ci) „2 f ^ 35 , A 



2a3sinh-'^(Ci) V 24/302^2^ 
where /(i) and J(t) time-dependent amplitudes are the following elliptic integrals, 

I{t) = 2-2/3V3A f '"'^"f^^""^ dr , J(t) = -^/(i) + 3sinh-i/^(Ci)cosh-i(Ct) , (15) 
Jo cosh (Or) 2^'-' 

which can be brought to the Legendre normal form given in the Appendix of 0]. The Sap represents gravitational 
waves analogous to the term containing Dap{x^ ,7]) in ([8]). We note that in Eqs. p2|) . ((T3)) and ((T5)) the integration 
functions Cq and J{t) are defined with an opposite sign compared to the original results of Perjes et al. in [6|. This 
sign convention helps us in the following to show the clear correspondence of the forthcoming results derived parallely 
from the Sachs & Wolfe and the A ^ solutions. 

In the chosen comoving coordinate system, the observers move together with the perturbations of the fluid in the 
perturbed spacetime. Using the expression for the four velocity given in ([7]) we can calculate the rotation tensor 
LUab defined by Ehlers Q 

UJab = ^f'J'liUc-d - Ud-c) , (16) 

where t°jj — (5°;, — u°'Ub and the semicolon denotes the covariant derivative with respect to the perturbed metric gab- 
The four velocity Ua is connected to the hoa components of the metric perturbation from which - with the use of the 
solution of the field equations in the case of vanishing cosmological constant - we have the following expression for 

Woo — , (^Oa = , UJaP — 7^V^(C/3^q — C'q,/?) • (17) 

Repeating the calculations in the case of A 7^ we have 

Woo = , Woq = , UJaP — -^V^(C/3,Q — Ca.p) ■ (18) 

The relations above demonstrate that the perturbed vorticity tensor is completely determined by the V^Cq. integration 
function. 



III. THE SACHS - WOLFE EFFECT FROM ROTATIONAL PERTURBATIONS 

In present section we calculate the contribution of rotational perturbations to the fluctuations of the cosmic mi- 
crowave background radiation. As a consequence of the Sachs - Wolfe effect [l] , the perturbations produce contribu- 
tions to the temperature fluctuations determined by the following relation 

with 

r] = rjR — w , = e^w , e^ep = — 1 , e'^ = const. (20) 

Here the vector e" represents the spatial direction of the light signal as seen by an observer moving together with the 
fluid and the E subscript denotes the instant of time of the emission. Applying the following identity for an arbitrary 
smooth function / = /(a;°) 

f,a-^dw = Ue"dw - fdw , (21) 
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FIG. 1: The time dependence of the (j}{t), ipit) and (^{t) functions in the presence and the absence of A. The subscript A denotes 
the corresponding functions in the model with cosmological constant. The vertical lines belong to the following events: tdec(A) 
and tdec denote the time of decoupling in the models with and without cosmological constant, while the fjj(A) and tn denote 
the present age of the universe in the two models respectively. 



and calculating only with rotational perturbations (terms containing Cq in the perturbed solutions), (jl9p can be 
written in the following general form, 

— = [(^(77) - ^(Ty) V^) ae"] I'' - mCo^e^dw , (22) 

were the (p, ip and ^ time dependent amplitudes are the following in the cases of vanishing and nonvanishing cosmo- 
logical constant respectively: 

48 2 1Q2 

m = - , ^iv) = - , av) = — , (23) 

rj^ T]^ T]^ 

: 8C2a2 , 1 ^ 32C^al cosh(Ct) 

3smh^'-^{Ct) smh^'''(Ci) 9 smh^' -^{Ct) 

In FIG.l we have plotted the (j), ip and ^ time dependent amplitudes on a logarithmic scale. The shape of the curves 
show that the presence of the cosmological constant results only about a ~ 1/3 factor difference from the single fluid 
model at each functions. This means that the fluctuations represented by rotational perturbations in the microwave 
background are ~ 3 times smaller in the presence of A. It is also remarkable that the magnitude of the amplitudes 
are decreasing exponentially since the time of decoupling, and the present day contributions are negligible comparing 
to early epochs. Except the integrated term in (j22p - whose contribution is growing in time ~ the contributions of the 
first two terms are the difference of their values at the time of decoupling and today. Therefore the dominant parts 
come from the time of decoupling and - as the present values are many orders of magnitude smaller - these terms 
produce the only relevant contributions to the temperature fluctuations that could be measured today. 

In the remainder of this section we present the recipe of the calculations how to obtain the anisotropy spectra of 
the fluctuations generated by rotational perturbations. As a first step, following Bardeen's covariant formalism, we 
replace the Ca{x^) perturbation function with a single mode of its decomposition into the vector harmonics Qa"'' of 
the Helmholtz equation [8| 
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where the subscript k denotes the wavenumber dependence of the coefficients C^'"'' and 

Q^°'> = exp(ik-x), 

Qi") = -fc-ig(2 , (26) 
Qi^^^ = -^(ei ± ie2)aexp(ik • x) . 

The spatial unit vectors ei and 62 are spanning the plane perpendicular to k. 

Applying this expansion, the integrated term in (j22p can be written in the following form in the case of vanishing 
cosmological constant 

isw = -^E^i^'^(-i±-2)"^"X ^3^^^" ■ 

Using the kfje^ = fccos(C) = n replacement one can evaluate this ISW integral in a closed analytic form as follows, 

ISW - - ^fc^'^ (^1 ± ^«'2)ae" [e™(i-"-' (6 - 2inr^E - n^E + i^^l) 

^^Ve ±1 

+ ryln^'e'" (r(0,m77£;)-r(0,m))- 7/1(6- 2m -n^+m^)] . (28) 
Repeating the calculations in the case of A 0, the integrated term in takes the following shape 

167(72^2,-^ fVE pin{l+x{y)) 

ISW (A) = i^^J^cf ^)(ei ±.e2).e" / dy , (29) 

where y = a/ao = sinh^^^^-* (Ci) and 



= ^_ r dy' _ /3 v/3-V3» ^/ /(3 + V3z)y 



SC'ao Jo + y'^ V 2 Cao 1 y 2(y + 1) 



(30) 



The ISW(A) integral in ([29]) can not be evaluated analytically thus numerical integration needs to be used. 

To be able to obtain the anisotropy power spectra generated by rotational perturbations, one has to evaluate the 
following well-known integral 

2i + i^^^vr-dKs\e^ , (31) 



47r 27r2 k 2i + 1 

where Ci is the anisotropy momenta and Qi{k,ri) is identical with 5T/T in (j22p - in the decomposition by vector 
harmonics of the Helmholtz equation presented above. 

As the last step, one has to take the square of the absolute value of 5T/T given in (HH). To be able to evaluate 
the integral expression of C^, one needs to have information about the spectrum of Ck which might be originated 
from predictions of models describing the universe before the decoupling era. As an example this has been modelled 
by Grishchuk in where he derives the angular variation of the microwave background radiation from rotational 
perturbations of quantum-mechanical origin. 



IV. ESTIMATIONS OF THE MAGNITUDE OF ROTATIONS 



As an application we estimate the magnitude of the angular velocity generated by rotational perturbations using the 
results of present observations of the CMBR fluctuations. It is well-known from the analysis of recent measurements 
Q that the dipole component of the temperature fluctuations has the largest amplitude with the order of magnitude 
~ 10^'^. This term certainly contains many contributions, however one of these originates from the effect of rotational 
perturbations. As we have shown in \22\ this contribution can be generally separated into the following three 
components: one proportional to Cq, one proportional to V^Cq, and an integrated term proportional also to Ca- 
Since the rotation tensor is completely determined by the V^Cq, function, we must estimate the order of magnitude 
of the V^Cq perturbation to be able to obtain the magnitude of the corresponding angular velocity vector Lj{t^x'^). 
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Applying the results of present observations of the dipole component, we can make an estimation on the upper bound 
of the rotations with the following procedure. 

The corresponding angular velocity deduced from the vorticity tensor is 

^ e'^^'ujps , (32) 



2^/7 

where 7 — — det(5Q/3) and e"'^^ is the 3-dimensional totally antisymmetric Levi-Civita symbol. Having done the 
calculation for the square root of cj^, we obtain the expression 

Lu{t, x") = w(t)K(x") with K(a;") ^ [(V^C^./j)^ - V^C^'^V^C/j,^] ^'^ , (33) 

and the w{t) time dependent amplitude is the following in the cases of vanishing and nonvanishing cosmological 
constant respectively, 

Hr 1 

"^^^ ^ V' "^^^^ ^ 2aosinh(V3)(^,) ■ (34) 

As an upper limit, in accordance with the measurements, we can say that none of the terms in (|22p can be larger 
than the order of magnitude lO""^, so this must also be valid for the V^Cq, term. In the previous section we have 
shown that the relevant contribution of this term comes from its value at the time of decoupling, where the magnitude 
of its time dependent amplitude (j){tdec) is known. This makes the restriction to W^Ca that its order of magnitude 
can not be larger than 10~^. If we assume that the first derivative V^Cq./j has the same order of magnitude, we 
can make the estimation for k{x°') in (j33p that is k{x°') 10~^. Applying this assumption we obtain the following 
order of magnitudes for the values of the angular velocity at the time of decoupling and today in the presence and 
the absence of A: 



uj^tdecx""), uJA{t dec, x'') ^ 10 "/year; uj{tB., x"'), ujA{tR, x"') ^ 10 /year . (35) 

It shows that the presence of the cosmological constant doesn't change the order of magnitude of the angular 
velocity, however in a more accurate calculation there might be a multiplying factor (< 10) between the two models 
similar to the results obtained for 5T /T in the previous section. 

The present value of w can be associated with the rotation of the universe, and the obtained order of magnitudes 
are in good agreement with the results of where the estimated value for the ratio uj{tii)/Hu < 10^^ in a flat 
model. 



V. CONCLUSIONS 



We have considered the effects of rotational perturbations of FRW models on the cosmic microwave background 
radiation and its connection to the amplitude of the rotation of the universe. We presented the general form of the 
relative temperature fluctuations induced by rotations and obtained analytic expressions for all the time dependent 
coefficients arising both in the presence and the absence of a cosmological constant. It is found that the and 
ip transfer functions of the rotational Sachs- Wolfe effect, and the ^ weighing function of the rotational integrated 
Sachs- Wolfe effect are decaying exponentially, and there is a ^ 1/3 factor between their values in the models with 
zero and nonzero cosmological constant. We presented the general method of calculating the anisotropy power spectra 
of the rotational ISW term, and using the results of recent observations for the dipole component of the microwave 
background, we made an upper estimation for the magnitude of the rotation of the universe at the time of decoupling 
and today. It is a remarkable result that the value of the corresponding angular velocity at the time of decoupling 
is close to the order of magnitude of the angular velocity of a typical galaxy cluster today [l^, [lll| . This interesting 
coincidence might deserve attention in future investigations. 
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APPENDIX: PHYSICAL PARAMETERS 

We present the list of physical parameters that have been used throughout the calculations obtaining numerical 
values. 



Hr = 70 km s-'^Mpc-^ ; (A.l) 

Q — VlMatter + f^A = 1 : ^Matter = 0.3 , = 0.7 ; (A. 2) 

A = SH^nA/c^ = 1.2 X lO^'^^m^^ ; (A.3) 

C = \/3A/2 = 2.85 X 10"^* s"^ ; (A.4) 

ao = 1.1 X 10^^ s ; (A.5) 

Zdec ^ 1100; (A.6) 

m = (A.7) 
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